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Abstract 

We use inelastic neutron scattering and molecular dynamics (MD) simulation to investigate the 
chemical short range order (CSRO), visible through prepeaks in the structure factors, and its 
relation to self diffusion in Al-Ni melts. As a function of composition at 1795 K Ni self diffusion 
coefficients from experiment and simulation exhibit a non-linear dependence with a pronounced 
increase on the Al-rich side. This comes along with a change in CSRO with increasing Al content 
that is related to a more dense packing of the atoms in Ni-rich Al-Ni systems. 
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In hard-sphere like metallic liquids mass transport is strongly connected to the packing 
fraction of the atoms . However, in liquid Al-Ni Al-Fe 0,0] and glass forming 

Al rich alloys [7] the structure clearly deviates from a random hard-sphere packing and 
exhibits a chemical short range order (CSRO) Q]. As was shown in a neutron diffraction 
experiment on liquid Al 8 oNi 20 |4|, CSRO manifests itself in a prepeak in the static partial 
structure factors at an intermediate lengthscale around 1.8 A -1 . In this paper, we study the 
composition dependence of the CSRO in Al-Ni melts and show how it affects the diffusion 
dynamics in this alloy. To address this issue, we performed inelastic neutron scattering 
experiments and molecular dynamics (MD) simulations. 

Al-Ni alloys were prepared by arc melting of pure elements under a purified Argon atmo- 
sphere. The corresponding solidus and liquidus temperatures were measured with differential 
scanning calorimetry and found to be in excellent agreement with the phase diagram 0. 
For the neutron time-of-flight experiment on the IN 6 spectrometer of the Institut Laue- 
Langevin, thin-walled AI2O3 containers were used that provide a hollow cylindrical sample 
geometry of 22 mm in diameter and a thickness of 1.2 mm for the Al rich samples, and 
0.6 mm for the AI25M75 sample. An incident neutron wavelength of A = 5.1 A yielded an 
energy resolution of 5E ~ 92 /ieV (FWHM) and an accessible wave number range at zero 
energy transfer of q = 0.4 — 2.0 A . Alioo-xNi^ alloys were measured above their liquidus 
temperatures at 1525 K and at 1795 K for x = 10,20,23,30, and at 1795 K for x = 38,75 
compositions. 

The scattering law S(q, u) was obtained by normalization to a vanadium standard, correc- 
tion for self absorption and container scattering, and interpolation to constant wave numbers 
q. Further, S(q, uS) was symmetrized with respect to the energy transfer fko by means of the 
detailed balance factor. Fourier transformation of S(q,u), deconvolution of the instrumen- 
tal resolution, and normalization with the value at t = gave the time correlation function 
<&(q,t). For times above ~ lps, <&(q,t) shows an exp(—q 2 D q t) decay to zero. Towards small 
q incoherent scattering on the Ni atoms dominates the signal and the g-dependent diffusivity 
D q becomes constant. Thus, Ni self diffusion coefficients D could be derived on an absolute 
scale 0,0. 

In order to model the interaction between the atoms in the MD simulation, we used a 
potential of the eroded ato m type that was r eeentl y derived by Mishin * a,. Q. The 
simulations were done at 1795 K and 1525 K for different Al-Ni compositions. The systems 
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consist of 1500 particles in each case. First standard Monte-Carlo (MC) simulations in the 



NpT ensemble 



12j were used to fully equilibrate the systems at zero pressure and to generate 



five independent configurations for MD simulations in the microcanonical ensemble. In the 
latter case, Newton's equations of motion were integrated with the velocity Verlet algorithm 
using a time step of 1.0 fs. From the MD trajectories the structure factors and diffusion 
constants were determined [13j. 

By integrating the scattering law S(q,u) over the quasielastic line inelastic neutron scat- 
tering also gives access to structural information on intermediate length scales. Fig. 1 
displays the integral over the quasielastic line (denoted by S qu (q) in the following) for the 
Al-Ni melts at 1795 K For clarity a constant that approximates the incoherent scattering on 
the Ni atoms has been subtracted. For Al-rich alloys the spectra exhibit a prepeak that is 
increasing with increasing Ni content and that shifts from ~ 1.6 A -1 in Al§oNi2o to ~ 1.8 A -1 
in Al62Ni3§. In AlgoNiio nearly no prepeak is visible anymore. As compared to the data at 
1795 K, the quasielastic structure factors measured at 1525 K do not reveal a measurable 
shift in the peak positions and the intensities in the maxima of the prepeaks exhibit a ~ 20 % 
increase which can be accounted for by the increased Debye- Waller factor. On the Ni rich 
side the Al 2 5Ni 75 S qn (q) also exhibits a prepeak. This prepeak is much broader, though, and 
has its maximum at ~ 1.5 A -1 . 

Fig. 2 shows the measured Ni self diffusion coefficients as a function of composition at 
1525 K and at 1795 K. At 1795 K values range from 3.95 ± 0.10 x lO^mV 1 in Al 25 Ni 75 
to 10.05 ±0.11 x lO^mV 1 in Al 90 Nii . The Ni self diffusion coefficient in Al 2 5Ni 75 is 
equal within error bars to the value in pure liquid Ni with D = 3.80 ± 0.06 x 10~ 9 m 2 s -1 
j3j. Substitution of 25% Ni by Al does not affect the diffusion coefficient significantly. In 
contrast, on the Al rich side diffusion coefficients show a pronounced increase with increasing 
Al content. 

The temperature and concentration dependence of the measured diffusion coefficients are 
well reproduced by the ones of the MD simulation: There is an overall 20% agreement with 
the experimental data (Fig. 2). From the inset of Fig. 2 we see also that in our simulation the 
ratio D^j/D^ does not depend on temperature in the considered temperature range (here 
the star means that diffusion constants from the simulation are considered). Note that the 
diffusion constants, that were obtained in a recent simulation study by Asta et al. ^| using 
different EAM potentials, are a factor 2-3 higher than those obtained in our simulation, but 
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also display a pronounced non-linear behavior as a function of composition. As we shall see 
in the following, the non-linear behavior of the diffusion coefficients is reflected in structural 
changes that can be understood from the detailed informatiom provided by the simulation. 

Fig. 3 displays the S n (q) = NAlti i l Nm ^ J2k,i b k k(e i(l - {rk ~ ri) ), with r k being the position 
of the fc'th particle. N a , a G Al, Ni, denotes the number of atoms of type a, and &ai = 
0.3449 • 10 -12 cm and = 1.03 • 10~ 12 cm are respectively the neutron scattering lengths of 
Al and Ni. As Fig. 3 shows, also a prepeak is observed in the simulation which is located 
at slightly smaller q for the Al rich systems whereas for AI30M70 the prepeak is at a slightly 
larger q than expected from the experiment. However, in particular for the Al rich systems 
a similar behavior as in the experiment is found, i.e., with increasing Al concentration, the 
prepeak becomes broader, its amplitude decreases and its location shifts to smaller q. We 
note that S qu (q) and S n (q) are different quantities, since S qn (q) is also affected by the q 
dependence of the Debye- Waller factor. Moreover, S qu (q) contains a mixture of incoherent 
and coherent contributions that cannot be disentangled from each other. Thus, one should 
compare only the relative intensities of the peaks in the experimental S qu (q). 

Also shown in Fig. 3 is the main peak which indicates repeating units of neighboring 
atoms. The location of this peak moves to higher q towards the Ni rich systems which is due 
;o the fact that with increasing Ni concentration the packing density of the atoms increases 
13. This can be inferred from the inset of Fie;. 3 which shows the atomic volume V a as 

n 

found in the simulation in comparison to experimental data [15J . Note that we get an overall 
5% agreement with the experiment over the whole Al concentration range. The nonlinear 
behavior of V a is strongly correlated with the nonlinear behavior of the diffusion constants 
clS cL result of the CSRO. 

S qn (q) and S n (q) are total structure factors which are in binary systems linear com- 
binations of three partial structure factors. In order to understand the CSRO that is 
present in Al-Ni melts one has to consider these partial structure factors defined by 
= y/NcNp ^£=1 Z^li( ex P(^l • ( r k a - r i /3 )))> where the indices k a and l p corre- 
spond to particles of type a and /3, respectively. In Fig. 4 the different S a p(q) are shown 
again at T = 1525 K for Al rich compositions, for a 50:50 mixture, and for the Ni rich com- 
position Al 3 oNi 70 . The emergence of a prepeak in S a p(q) indicates that there are repeating 
structural units involving next-nearest a(3 neighbors which are built in inhomogeneously 
into the structure. As we see in Fig. 4, no prepeak emerges in Saiai(<z) f° r f ne Al rich 
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compositions but the main peak moves slightly to higher q with decreasing Al concentra- 
tion. By further decreasing the Al concentration a prepeak evolves that moves to smaller 
q. As a result a double peak structure is observed in Saiai(<7) f° r Al 30 Ni 7 o with maxima at 
qi = 1.8 A -1 and q 2 = 3.0 A -1 . In Sj^n(q) a prepeak with a negative amplitude is visible at 
all Al concentration and, as in the Al-Al correlations, this prepeak moves to smaller q with 
decreasing Al concentration. The negative amplitude of the prepeak in SainK?) is due to 
the avoidance of corresponding Al-Ni distances. Also in Smm(q) a prepeak is present at all 
concentrations. This prepeak broadens with increasing Al concentration while its location 
remains essentially constant. In AI30NI70 the main peak is at the same location in the three 
different partial structure factors, namely around q = 3.0 A -1 . 

The picture provided by the simulation is as follows: In the Ni rich systems the minority 
species Al is built into the structure such that the distance between repeated Al-Al units 
corresponds to that of repeated Ni— Ni units [16] whereas towards the Al rich systems the 
distance between repeated Al-Al units becomes larger without a significant change in the 
distance between repeated Ni-Ni units. All this explains the nonlinear behavior of the 
atomic volume and, since it indicates a nonlinear behavior of the packing of the atoms upon 
a change of the composition, it demonstrates that the nonlinear increase of the diffusion 
constants upon increasing Al content is essentially a packing effect. Of course, this packing 
argument does not explain the deviation of D* M / from one for the Ni rich systems. The 
understanding of this feature requires a microscopic theory. 

In conclusion, an investigation of Al-Ni melts using inelastic neutron scattering and 
MD simulation has shed light onto the nonlinear composition dependence of the diffusion 
constants and its relation to chemical short range order visible through a prepeak in the 
structure factor. We confirm the finding for other hard-sphere like metallic melts [3| that the 
transport in these systems is strongly affected by packing effects although the interactions 
between the atoms depend strongly on chemical details. 

We are grateful to Kurt Binder for stimulating discussion. The authors acknowledge 
financial support by the Deutsche Forschungsgemeinschaft (SPP Phasenumwandlungen in 
mehrkomponentigen Schmelzen) under Grant No. Bi314/18 and Mel958/2. One of the 
authors (J. H.) acknowledges the support of the DFG under Grant No. HO 2231/2. 
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FIG. 1: Quasielastic coherent structure factor S qu (q) at intermediate q of Al-Ni melts at 
1795 K as obtained by inelastic neutron scattering (see text). 

FIG. 2: Diffusion coefficients at 1795 K and at 1525 K as a function of the aluminium 
concentration from simulation and neutron scattering. The broken lines are guides to the 
eye. Inset: Ratio of Al and Ni self diffusion coefficients from simulation. Note that the 
diffusion constants from the simulation are decorated by a star. 



FIG. 3: Structure factors S n (q) of Al-Ni melts as calculated from the partial structure 
factors of the simulation weighted with the neutron scattering length. Inset: Atomic 
volume from the simulation compared to experimental data from Ref. 3] as a function of 
composition. 



FIG. 4: Partial structure factors Sa^q) of Al-Ni melts as determined by the simulation at 
1525 K. a) Saiai(?), b) S A im{q), and c) Smm(q)- 
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